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C^ ' Abstract 

o , 

QQ We discuss, in the framework of perturbative QCD at next to leading order, two related 

f^ , observables which are usually considered to provide tests of the BFKL dynamics : jet-jet 

correlations at Tevatron energies and forward particle-jet correlations at HERA. In the first 

case we study the rapidity gap dependence of the azimuthal correlations and find slightly too 

^^ ■ strong correlations at large gap. In the second case we discuss the cross section as well as the 

$_( ' azimuthal correlations over a rapidity gap range of 5 units. We find that the requirement of a 

5t 1 forward particle imposes strong kinematical constraints which distort the distributions, notably 

at small rapidity gaps. We also show that the decorrelation is stronger in electroproduction 
than in hadron-hadron collisions. Unfortunately no data are yet available for comparison. 
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1 Introduction 

The first test of BFKL IJ dynamics in liadronic collisions was proposed by Mueller and 
Navelet p] who showed that the cross section for the inclusive production of 2 jets at the partonic 
level would increase exponentially with the rapidity separation, A77, of the jets. At the hadronic 
level the cross section involves a convolution with the initial parton distributions, as a result of 
which the hadronic cross section decreases with the rapidity gap due to the very fast fall-off of 
the parton densities. In order to disentangle the interesting dynamical effects from the parton 
distributions effect the study of the normalized azimuthal angle distribution [31 [3] was proposed 
which should be less sensitive to the latter. More precisely, the relevant quantity to consider is the 
average < cos(7r — 0) >, where (j) is the azimuthal angle between the 2 jets. From BFKL dynamics 
one would expect a very fast decorrelation, due to the emission of a growing number of jets in the 
rapidity interval between the observed jets. This is to be contrasted to the case of DGLAP [5] 
dynamics as implemented in lowest order (LO) and next-to-lowest order (NLO) calculations: at 
LO the 2 jets are perfectly correlated (0 — tt) while at NLO the decorrelation arises from the 
emission of just one extra jet. Monte-Carlo codes based on DGLAP dynamics offer however the 
possibility of a stronger decorrelation due to parton showers. 

Coming back to BFKL dynamics it was noted [6l[71[8] that the predictions were very sensitive to 
energy momentum constraints neglected in the original analytic calculations. This is, in particular, 
true for azimuthal correlations where an exact energy momentum conservation considerably softens 
the rapidity dependence, compared to the analytical (asymptotic) results. Indeed, less phase space 
is available for the emission of intermediate gluons in the former case. Because of these and 
other effects [U |9] the BFKL predictions show less decorrelation than that found in the initial 
calculations. 

Data on < cos(7r ~ cj)) > a.t Tevatron energies have been obtained by the D0 collaboration [TD] 
in the range < Arj < 5 and have been compared to a variety of calculations. The leading order 
BFKL asymptotic predictions are completely ruled out by the data while the leading order ones, 
including energy momentum conservation [SJ [5] , or the next-to-leading order ones [5] are closer 



to the data with still a too strong rapidity decrease. On the other hand, the NLO predictions of 
JETRAD [Tl] do not show enough decorrelation as the rapidity increases. Only the HERWIG [12] 
results are in very good agreement with the data. 

BFKL dynamics can also be probed by considering the cross section for forward production of a 
large p^ jet or particle in deep- inelastic scattering (DIS) events at small Bjorken-a; xbj il3j, in the 
domain p^ ^ Q^. Such data have been collected by the HI [HI [15] and ZEUS [16: collaborations. 
For this observable, models based on BFKL dynamics (supplemented by some phenomenological 
inputs) , 17] , were claimed to be in reasonable agreement with data although they did not account 
for the correct Q^ evolution. More recently, several NLO calculations of forward 7r° production 
have been carried out [18 ] fT9 ] [20 ] [211 [22] [23] and found in agreement with the latest HI data [15] : 
in particular, the dependence of the cross section as a function of the variables xbj, Q^, x-,^,p*^ 
(where a;^ and p*^ are, respectively, the fraction of longitudinal momentum and the transverse 
momentum of the pion in the 7*-proton center of mass frame) are very well reproduced [20J. 

To further study the validity of the perturbative approach in deep-inelastic scattering, we pro- 
pose to consider azimuthal correlations involving a forward large p±^ hadron, to probe the small 
XBj where BFKL dynamics is often thought to be more relevant. Because of the different subpro- 
cess structure, it is interesting to compare the azimuthal correlations in hadroproduction and in 
electroproduction. However, in the latter case, the perturbative dynamics is more complex because 
the cross section is built from two contributions, the so-called direct and resolved contributions. 

In the following section, we first go back to correlations in p p collisions at the Tevatron, from a 
perturbative point of view, and find that the rapidity dependence, quite sensitive to scale changes, 
is within the experimental errors, although being not steep enough. In Sect. [3] we turn to DIS 
and study forward-pion jet correlations. We consider two classes of observables: the cross section 
and the average < cos(7r — (/))> as a function of the rapidity hadron-jet separation for different 
Q^ values. The importance of kinematical constraints on the shape of the rapidity distribution is 
stressed. Conclusions are given in a final section. 



2 Azimuthal correlations in hadron-hadron collisions 

The D0 collaboration measured [Tp] the rapidity dependence of the azimuthal correlations 
between large-pj_ jets at ^/s — 1800 GeV with the following criteria: among the jets with E±_ > 
20 GeV, the two jets having the largest and the smallest pseudorapidity were selected. One of 
the two jets was also required to be above 50 GeV. The selected jets are required to be in the 
pseudorapidity range |ryi,2| < 3 and the azimuthal correlations are measured as a function of 
A77 ~ Tji — ri2- More precisely, D0 measures < cos(7r — (/))> as a function of A77 where cf) is the 
azimuthal angle between the two jets. 
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Figure 1: Theoretical predictions, at NLO in perturbative QCD, of the two-jet correlations 
(squares: scale P±/2\ diamonds: scale i'±/4) compared with D0 data |10: (solid dots). The 
solid line indicates the systematic errors. 



To study this quantity we use the code DUET which belongs to the PHOX Family [SJ and 
calculate single jet and dijet cross sections at NLO in hadron-hadron collisions. The interested 
reader may find technical details on the PHOX Family codes on the site 24 . We have checked 
that the DIJET predictions are in good agreement with the single jet cross section of ref. [33] 
and with the spectrum dN/dA-q measured by D0 [10] under the conditions described above. The 



parton distribution functions are CTEQ6M [26] and we use Pi/2 as the standard factorization 
and renormalization scales, where P± is half the sum of the transverse momenta of the largest p± 
partons, Pi = £J=iiHJ=a. For the jets we use Re = .7 and Rsep = 1-3 Re as in ref. [25] , 



We now focus on the azimuthal correlations. We calculate the values of < cos(7r — (j>) > quoted 
by the D0 collaboration in Fig. 4 of ref. [TOj. We expect the corresponding theoretical predictions 
to be sensitive to the choice of scales. Indeed the aplanarity of the jets is due to the 2^3 
subprocesses and when is different from tt we perform a LO calculation. Only the region (p ^ tt 
corresponds to a full NLO calculation. To test this sensitivity we start with the standard scale 
P±/2. The result is displayed in Fig. [1] (squares) together with the D0 data (solid dots). The 
theoretical error due to the limited statistics in the Monte Carlo integration is of the order of 
± 10 % or less. 

Another scale which may appear more natural is (minus) the square of the momentum transfer 
in the t-channel of a 2 — > 2 subprocess, Q^ — p'^^{l + e^*^''^^''^)), pj^ being a jet transverse momen- 
tum. In order to have a scale close to p±/2 when rji ^ 772 we use c P± ( ^-^ — g" — ~ ) with c = i 
for the factorization and the renormalisation scales. We obtain results almost identical to those 
obtained with the scale P±/2 (about 1 % lower). Finally with the scale Pj_/4 we get the results 
shown by the diamonds in Fig. [T] 

The agreement between data and theory is not perfect, the theoretical points staying at the 
upper limit of the experimental error (statistical + systematic). The smooth decrease of < cos(7r — 
(f)) > with increasing A77 is well reproduced at small and medium rapidity gaps, but at large A77 the 
jets in the NLO calculation appear to be too much correlated. Note that with the scale X)i=i P-Li/2 
we recover the JETRAD results quoted in [TD]. 

Other approaches based on the generator Herwig [12] and on a BFKL calculation [6] are also 
displayed by the D0 collaboration. The Herwig results are in good agreement with data whereas 
the BFKL prediction undershoots them by a large amount. A recent BFKL analytic calculation 
at next to leading order |9] gives a somewhat better description than the leading order one of Del 
Duca and Schmidt quoted in lOj which takes into account some non asymptotic effects. 

The conclusion of this study is that HERWIG, and to a lesser extent NLO calculations, can give 



a satisfactory description of the data, although the underlying dynamics is different. In the NLO 
calculation the jet aplanarity comes from the 2-^-3 subprocesses described by exact QCD matrix 
elements. These matrix elements are not present in HERWIG where the aplanarity comes from 
the initial and final parton showers. This mechanism leads to a correct value of < cos(7r — 0) >. 
(However the cross section for three large-pj^ jets is expected to fall an order of magnitude below 
data, because of the lack of exact 2 — > 3 matrix elements in this approach). Based on this analysis, 
we speculate that a NNLO calculation would give a satisfactory description of the data. 

3 Hadron-jet azimuthal correlations in DIS 

In ref. [20j we carried out a detailed study of single inclusive forward tt'^ production in deep 
inelastic scattering and compared the NLO perturbative results with the recent HI data (see 
ref. [m [15] for experimental details and cuts). The data probe a small xbj range, 1.1 10^* < 
XBj < 1.1 10~^ where DGLAP dynamics may break down. The calculation includes the "direct" 
component (where the virtual photon couples to the hard processes) and the "resolved" one where 
the photon acts as a composite object, both contributions being calculated at the next to leading 
order. We use the FORTRAN code DISPHOX developed from the real photoproduction code 
JETPHOX :24] and described in refs. [HI [TH [20] with the following inputs: the proton structure 
functions are taken from the CTEQ6M tables [26], the parton fragmentation functions into pions 
are those of KKP [27] while, for the resolved case, the photon structure functions are found in [TO] . 
We work in the MS renormalization and factorization schemes and all the scales are equal to 
(Q^ +pi). We take n/ = 4 flavors and A^g^ = 326 MeV. 

In Fig. [5] we display the single inclusive forward pion distribution as a function of xbj and its 
comparison with the various theoretical components of the cross section. A remarkable fact is the 
smallness of the Born direct term and the very large size of the higher order correction to this 
Born term (see [20 for a precise definition of various components). As already noted, at lowest 
order, the hard processes are mediated by quark exchange while, at higher orders, there appear 
diagrams with a gluon exchange in the t channel, e. g. 7(7 -^ qqq: because of the gluon pole these 
terms are enhanced in the case of forward production of hadrons and they considerably increase 
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Figure 2: The cross section da/dxBj corresponding to the range 4.5 GeV^ < Q^ < 15 GeV^ 
compared to HI data [I^. The various theoretical components of the cross section are shown. The 
symbol HOs denotes the HO correction from which the lowest order resolved contribution has been 
subtracted. 





Figure 3: Examples of higher order corrections to (a) the direct contribution, (b) the resolved 
contribution. The forward hadron is emitted from the bottom parton line. 



the production rate. The importance of the resolved component can also be explained along the 
same lines because, already at lowest order, it contains processes with a gluon exchange. In that 
case, the hard processes, at LO and NLO, are the same as in pp collisions, the only difference with 
the latter case being the convolution with the structure functions: those of the photon are less 
steep than those of the proton and the relative weight of gluon and quark initiated processes are 
different. In Fig.[3l we show two diagrams characteristic of the higher order corrections included in 
our calculation: it is clear that these diagrams can be seen as the lowest order terms of the BFKL 
ladder stretched between the photon (backward direction) and the proton (forward direction). If 
the available phase space (rapidity) between the forward and backward partons is large enough 
then multiple gluon emission is expected and the BFKL dynamics should dominate. In the opposite 
case the emission of one (or two) gluons is probably more appropriate. The very good agreement 
between the inclusive HI data , l5j and our perturbative calculation |20j seems to support the latter 
hypothesis which can be further tested by considering forward hadron-jet azimuthal correlations. 
The HI collaboration has studied azimuthal correlations between the highest two transverse 
momentum jets |28] in DIS events at small xbj ■ They find a reasonable agreement with perturbative 
calculations except for small separations in azimuth where models based on BFKL dynamics fare 



better. Unlike the analysis in p p collisions they do not study the correlations between largest 
rapidity gap jets nor do they measure < cos(7r — </))> as a function of the rapidity gap. 

In the following we study particle-jet correlations and we use the same DISPHOX code which 
was successfully used to describe the single inclusive forward pion data and we concentrate on 
the rapidity dependence of cross sections and < cos(7r ~ (f>) >. Electron-jet correlations have 
already been considered in [291 ISQ] in the framework of the BFKL approach. In contrast to p9l 
I30j . we propose to measure the forward hadron-jet correlations in which the trivial effect of the 
electromagnetic vertex is not present. These correlations are also closer to those measured in 
hadron-hadron collisions. 

The momentum of a large-pj^ final state is better fixed in experiments measuring hadrons than 
in those measuring jets. In the latter case, the link between jet and parton momenta are obtained 
via Monte-Carlo calculations and it might not be well controlled for small transverse momentum 
jets, with p± < 4 GeV/c say. Besides, using a particle rather than a jet also allows us to extend 
the useful rapidity range of the observable. 

All our results will be presented in the hadronic center of mass frame (CM*) where the virtual 
photon and initial parton in the proton are aligned along the z axis. The reason is the following. 
In the laboratory frame the photon has a transverse momentum {qf = (1 — y)Q^X!j and therefore it 
can emit a collinear parton which can be part of a large-pj^ jet. In the resolved case this parton is 
also collinear to the spectator partons contained in the photon and these latter must also be part of 
a large-p_L jet. As a result, large-p_L jets defined in the laboratory frame can be made of spectator 
partons, and we would need to modify the jet algorithm used in DISPHOX which only considers 
jets made of active partons (produced in the subprocess). Besides, the cross section would be 
very sensitive to the position of the cone (defining the jet) with respect to the photon momentum. 
For instance, a cone tangent to the photon momentum will lead to a "collinear divergent" (the 
quotation marks indicate that the divergence is actually regularized by the photon virtuality) 
cross section and factorization is broken. In order to avoid such pathologies, we never work in the 
laboratory frame and we impose, to insure the applicability of perturbation theory, a minimum 



tj; is the inelasticity of the electron y = q ■ P/Pe ■ P, with Pe and P the initial electron and proton momenta 
respectively. 



value on the transverse momentum of particles measured in the CM* frame. 

For our calculation we use the following cuts, in agreement with the HI detector constraints. 
The lepton inelasticity is limited by 0.1 < y < 0.7 and the Bjorken variable by 1. < IC* xbj < 
5. The forward tt" is observed in the laboratory in the domain 5° < 6^^^ < 25° with a scaled 
energy .01 < x^ = E^^^/E^p^ < 1. Furthermore we impose, in the CM* frame, p*^"^ > 4 GcV and 
p-^jct > 3 5 QeV. This last cut is defined to avoid too large higher-order corrections in some regions 
of phase space. The range of Q^ is limited to 50 GeV^ and, in our study, we vary the lower bound 
from 2 GeV^ upward. 

Before presenting the results it is worthwhile to discuss the constraints imposed by the above 
cuts. First, let us recall the necessary conditions for the BFKL dynamics to be relevant: one should 
maximize the rapidity extent of the photon-parton system. If W be the energy of this system, we 
should maximize the quantity (we recall the relation yS — Q^/xbj)- 

W^ _ \fW (a^i - ^Bj) _ VyS {xi - XBj) 



(3.1) 



where S is the total electron-proton energy squared, xi the fraction of the proton momentum 
carried by the struck parton and p*^ the transverse momentum of the final state parton. As is well 
known, this is achieved by using large y and small xbj, for p*^ not too large. Considering the Born 
term kinematics with two hard partons in the final state, we can derive the following relation in 
the CM* frame 



Q^(l - xi) _ xi - xbj 
yS - Q"^ 1 - XBj 



exp(2r/^ - A77 ) = xi - —-^ — -5- = -; ^ ~ xi - XBj , (3.2) 



where 77* is the rapidity of the observed forward pion in the CM* frame and A77* the rapidity 
gap between the pion and the jet. The angular conditions on the forward pion constrain the pion 
rapidity in the laboratory, 1.50 < 77^^'' < 3.13, or approximateljo in the CM* frame —.65 < 77^ < 1. 
When looking at the pion-jet correlations, if we further impose the condition of a forward going jet, 
Ary* decreases and xi increases. Therefore the condition A77* small corresponds to a forward jet 
implying a large value of Xi where the proton structure function is suppressed. Another relation 



tin the CM* frame, the pion rapidity is approximately (Q^ ■C P*^) given by r;* = r]^^^ — 1.7+ ^ In 1/ and we have 
used i In J/ ^ —.46 for j; ~ .4 



is (for XBj <C xi) 



„2 „„„i,2 ^V* _ Q^ (a:i -xbj) _ XiyS -Q"^ S 



Pl cosh' -f- . t ""^^ - ^^"l-^ ~ -" 4 ^ <"' 



showing that at large values of A77* correspond large values of y (and small values of xi by eq. 13. 2p 
which are cut by the condition y < .7. 

These features can be clearly seen in Figs. |4] where the jets are defined as those with the 
largest transverse momentum. No rapidity plateau appears in the A77* dependence of the cross 
section and the maximum of the cross section moves to larger values of A77* when Qmm^ ^^'^ 
correspondingly y, increases, as expected. In these figures it is worth noting that the resolved cross 
section, the relative weight of which decreases with increasing Q^ , is less affected by the kinematics 
constraints discussed above in the Born case. Because it involves an extra convolution with the 
parton distribution functions in the photon, it is flatter than the direct contribution. 

Some of the kinematical issues can also be understood by writing y and xi in terms of the 
rapidity of the final partons j/i and y2- In a frame in which the proton and photon are coUinear 
(P is the proton momentum in this frame) 

n* IP 

and 

x, = ^^{ey^+ey^)+XB, (3.5) 

These equations directly give the variations of xi and y in terms of the rapidity of the final state 
jets and thereby the dependence of the parton distributions of the proton and the photon on these 
rapidities. Thus, a second forward parton corresponds to a larger value of xi whereas a backward 
parton corresponds to a larger value of y. These kinematical domains are thus suppressed by 
the large xi and large y behaviour of the proton and photon distribution functions respectively. 
Moreover, in the resolved case, ea. p.4[) is modified by multiplying the y on the LHS by an extra 
convolution z due to the further splitting of the photon into a qq pair. Thus, for the direct case, 
j/2 is determined by Q^, y and j/i which also determines the cross section. In the resolved case. 
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Figure 4: Variation of the cross section as a function of the rapidity gap between the forward pion 
and the largest p± jet. 



11 



because of the extra convolution as explained above, j/2 is shifted to the right, and hence the whole 
resolved cross section, after averaging over Q^, y, yi and z is also shifted to the right. 

There is a marked difference between the rapidity dependence of these cross sections and that of 
the dijet measured in ref. [28] . There, the dijet cross section is shown as a function of the rapidity 
gap between the two largest p± jets measured in the CM* frame but no special requirement is 
imposed on the rapidity of one of the jets. The cross section is maximum at A?;* = and 
slowly decreases with increasing rapidity gap. The difference in shape with our results is due 
to the requirement that the n be observed in the forward direction, in the laboratory frame, 
which introduces strong kinematical constraints which distort, in our case, the rapidity dependence 
of the cross section. This is confirmed by later studies of the HI Collaboration ^31] where the 
rapidity dependence of the cross section between a forward jet and a di-jet system of small rapidity 
separation, is shown. A peak is observed at A77* ~ 1 with the cross section rapidly falling on either 
side. Let us note that this observable does not allow us to cover a large rapidity gap domain (about 
2 units for the jet-dijet correlation instead of 5 in the particle-jet case). 

We now present in Figs. [S] the cross section da/dArj* where A77* is the rapidity gap, still 
measured in the CM* frame, between the forward pion and the jet (with a minimum transverse 
momentum as specified above) which has the largest rapidity gap with the pion. The same general 
features as in Figs. [Hare observed with however a somewhat fiatter behavior than before at medium 
rapidity gaps and a slower decrease at large gaps. This is related to the fact that the most backward 
jet is not necessarily the largest p± one. This observable can be compared to the one discussed by 
the D0 collaboration where two jets with the largest separation are considered. 

Let us now discuss some features of the NLO calculation. The HO corrections to the Born cross 
sections can be important in some kinematical domains. Consider, as an example, the "largest- 
rapidity-gap" case and the value (3^j„ = 10 Gev^. For the direct contribution the ratio HO/Born 
is equal to 1/3 in the range Arj* « 0-1 while, in the range A?/* — 4-5, it reaches 20. This large 
variation is due to the opening of a new channel, namely 7*3 —^ qqg (cf. Fig.[3ja)), which becomes 
important at large values of A77*. The variation of the ratio HO/Born is weaker in the resolved 
case: it varies from 1/3 to 4 for A77* increasing from the range 0-1 to the range 4-5. This increase 
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Figure 5: Variation of the cross section as a function of the rapidity gap between the forward pion 
and the jet with the largest rapidity gap with the pion. 
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Figure 6: Average < cos(7r — </))> as a function of the rapidity gap between the forward pion and 
the most backward jet. The different curves correspond to the minimum values of Q^ used in the 
calculation. 

is due to the larger phase space available for the third jet when Arj* increases. This behavior of 
the ratio HO/Born partly explains those of the average < cos(7r — </))> which we now consider. 

Figs. [SI QU give the dependence of the average < cos(7r — 0) > as a function of the rapidity 
gap under the same conditions as those of Fig. [5l The decrease of < cos(7r — </))> when A?7* 
increases is explained by the increase of the phase space available for the final jets in the 2^3 
subprocesses. Going into the details of the underlying mechanisms we see in Fig. [7] that the direct 
component shows a more rapid decorrelation as Ar/* increases because of the rapid growth of the 
HO contribution to this process. Comparing with the D0 data of Fig. [U we note that the overall 
size of < cos(7r —(/))> is smaller for the same value of Arj* but that the slope is rather similar at 
large Arj*. 

Coming back to Fig. [6] it is interesting to note the decrease of < cos(7r — 0) > when Q^j„ 
increases. This corresponds to the 2 — > 3 subprocesses becoming more isotropic. Finally the dip 



^Note that the errors in Fig. 7 are large primarily because of the errors in the resolved case. As usual the errors 
in the cross section and the cos </> weighted cross section have been added in quadratures to obtain the errors in the 
< cos(7r — (/>)> calculation. 
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Figure 7: Detail of the components of < cos(7r — 0) > as a function of the rapidity gap between 
the forward pion and the most backward jet for Q^ > 2 GeV^. 

at Arj* close to zero can be explained in the following way. First we note that in this A77* range, 
the HO corrections are smaller than the Born term. In the direct case, they even become negative 
when Q^in increases. In fact, this means that there are large compensations between the negative 
quasi 2 — > 2 HO contributions (virtual and collinear) and the positive contributions coming from 
the 2^3 processes. When the 2^3 contributions are weighted by cos(7r — (p), the (weighted) 
HO contributions strongly decreases, which causes a decrease of < cos(7r — (/))>. For higher values 
of Arj* , this effect is less pronounced. Let us also note that the resolved component is dominant in 
this bin for Q^ > 10 GeV^. Without HO corrections to this component, we would have obtained 
< cos(7r - (j)) >^ 1.0 

The behavior of < cos(7r — (/))> with A77* is strongly dependent on kinematical constraints, 
which are exactly respected in our NLO calculation, and this effect is specially strong near A77* ^ 0. 
On the other hand the BFKL approach docs not usually implement total energy momentum 
conservatiorlj and we expect quite different predictions. 

^For a discussion see ref. [8]. 
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4 Conclusion 

In this work we compare the azimuthal decorrelation effects, as a function of the rapidity gap 
between two jets at the Tevatron and between a forward tt*^ and a jet at HERA, using a next to 
leading order formahsm. In the case of hadronic colhsions, the calculated decorrelation effect is 
slightly slower than the data but it is still compatible with them taking into account statistical and 
systematic errors. In the deep inelastic case, the decorrelation is somewhat stronger than in the p p 
case, increasing with Q^, and we note interesting kinematical effects, related to the requirement of 
a forward particle trigger. In particular, it introduces a dip in the cross section at small rapidity 
gap, between the hadron and the jet. Dynamical effects and kinematical constraints also lead to 
very specific variations of < cos(7r ~ cj)) > with Arj*. 
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